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ATTRACT 


The  effect  of  chemical  kinetics  on  rocket  nozzle  design  was  in¬ 
vestigated  for  the  particular  case  of  hydrogen  gas  flowing  adiabati- 
cally  through  a  typical  rocket  nozzle  having  a  chanber-to-throat-urea 
ratio  of  2  to  1.  For  coi  parative  purposes  the  following  types  of 
flow  were  considered:  (a)  constant  composition  (frozen  equilibrium), 

(b)  instantaneous  chemical  equilibrium.  (shifting  equilibrium),  and 

(c)  kinetic  che:idcal  e  iuilibriun. .  A  stepwise  iteration  process  was 
employed  to  oerforn  the  integration  of  t:.e  differential  equation  in 
case  (c).  In  each  case  th.e  gas  entered  the  nozzle  at  a  temperature 

of  3500°K  and  a  pressure  of  20  atm  and  was  allowed  to  expand  isotropi¬ 
cally  to  an  exhaust  pressure  of  1  atm.  The  results  for  a  mass  flow 
rate  of  1000  grr.  sec  ire  presented  in  Figs.  2  through  6.  The  instan¬ 
taneous  equilibri  flow  assumption  not  only  gives  a  higher  specific 
impulse,  but  list  retires  a  larger  nozzle  than  either  the  Kinetic  or 
the  constar  position  flow  assumption.  I  he  kinetic  equilibrium 

flow  results  ire  intermediate  between  t! ose  for  instantaneous  equili¬ 
brium  flow  and  those  for  constant  composition  flow,  the  relative  posi¬ 
tion  depending  on  the  magnitude  of  t*  -  reaction  rate  whici.  governs  t  - 


INTRODUCTION 


It  has  been  standard  practice  in  the  firld  of  rocketry  to  compute 
the  specific  ir pulse  (i.e.,  the  jet  velocity  divided  by  the  accelera¬ 
tion  of  gravity)  of  a  propellant  system  on  tne  basis  of  either  of  two 
fundamental  assumptions,  namely,  that  the  .adiabatic  flow  of  propellant 
gases  from,  the  combustion  chamber  through  the  nozzle  is  such  that  either 
(a)  constant  composition  or  (b)  instantaneous  chemical  equilibrium  is 
maintained.  For  .an  adiabatic  process,  assumption  (a),  which  says  that 
the  inherent  reaction  rate  is  zero,  gives  a  mini mum  value  for  the  speci¬ 
fic  impulse,  vh-r-as  assumption  (b),  w.  icl  says  that  tie  inherent  reac¬ 
tion  rate  is  infinite,  gives  a  maximum  value  for  the  specific  impulse. 

The  introduction  of  reaction  rate  equations  into  the  computation  of 
specific  impulse  so  corr.pl i cates  the  problem  that  heretofore  investiga¬ 
tors  have,  perforce,  been  generally  satisfied  with  a  knowledge  of  the 
limiting  values  obtained  or.  the  basis  or  the  aoove  assumpt.i 

Aside  from  t!  •  computational  difficulties  w;  i ci  it  entails,  the 
kinetic  approacl  ,  because  it  requires  close  dr  **re-;ce  to  physical  reali¬ 
ty,  tends  not  only  to  establish  Ue  tru*-  specific  i --poise  for  a  propel¬ 
lant  c  stem.,  !  ut  tlso  to  give  some  insig''’  i  r.to  '  1,»»  d -sign  of  the  nozzle 
for  that  system.  Ir,  general,  specific  i  .pulse  value  cc  nut-d  ^nder 
ibsu  ptions  (aj  and  (b)  are  r.ade  wlti  out.  regard  to  u  zl*  configuration 
and  with  tie  added  simplifying  asriv  ptio?  th--t  Uir  chamber  aces  ^nter 
tie  nozzle  with  zero  velocity.  Diner  t : .  i s  lutt-r  arsu-ption  has  rr  ysb 
col  justification  .a:  d  is  erroneous,  it  leads  to  s  rious  v 


specific  impulse.  In  order  to  obviate  this  condition  it  is  nec-ssary 

at  t!  e  outset  to  determine  the  chamber  velocity  at  t.he  entrance  to  the 

nozzle  for  the  type  of  gas  flow  under  con.  .[deration.  This  may  be  done 

only  by  consideration  of  no/.zle  configuration. 

The  nozzle  generally  accented  by  rocket  designers  is  the  convergent 

divergent  de  Laval  type,  whicn  is  formed  from  two  intersecting  right  cir 

cular  cones.  The  half  angle  of  the  convergent  cone  is  usually  )0°  and 

that  of  the  divergent  cone,  15°.  To  prevent  discontinuity  of  flow,  the 

throat  section  is  rounded  off  with  a  r  dius  of  curvature  which  is  e  jual 

to  the  throat  diameter  Pr^.  Such  a  nozzle  was  conside-ed  by  airman  and 

Fenner.  The  nozzle  considered  in  this  study  is  shown  in  Fig.  1. 

It  differs  from  tie  41  tr  .an  an  i  Pinner  r.czzle  in  two  respects.  First, 

according.;  to  current  design  practice,  the  entrance  area  A  is  twice  the 

c 

thrc*  t  area  A  ,  This  restriction  peri. .its  one  to  establish,  among  other 

w 

things,  the  origin  of  a  rectangular  coordinate  system  at  the  entrance 

to  the  nozzle,  secondly,  the  convergent  cone  is  replaced  with  a  ,spi  ere 

of  radius  r  .  The  sphere  is  tangent  to  tie  rounded  throat  section  at  a 
c 

point  where  the  slope  rakes  an  angle  of  23°  31 1  with  the  nozzle  axis. 

Tie  use  of  a  spherical  section  instead  of  a  convergent  cone  ensures  con¬ 
tinuity  of  flow  from  th ■«  combustion  chamber  into  the  nozzle  entrance. 

The  propellant  system  in  this  study  is  pure  hydrogen  having  .a  cham¬ 
ber  temperature  of  3500°K  and  a  chamber  pressure  of  PO  at:;..  Its  history 
prior  to  its  arrival  at  t:.e  entrance  to  tie  nozzle  is  irrelevant.  From 
the  basic  nozzle  design  parameters  and  mass  flow  rat-  considerations, 
the  following  properties  may  be  computed: 


1.  Chamber  velocity  vc 

2.  Nozzle  throat  area  A^ 

3.  Nozzle  length  L 

4.  Jet  velocity  v^  or  equivalent  specific  impulse  I 

In  this  study,  the  above  properties  are  determined  not  only  for  (aj 
constant  composition  flow  and  (b)  instantaneous  chemical  equilibrium 
flow,  but  also  for  (c)  kinetic  chemical  equil Lbrium  flow.  Case  (cj 
r.  ikes  use  of  ti  e  expression  for  the  rate  of  recombination  of  hydrogen 
atoms  during  t:.e  adiabatic  expansion  process  of  the  propellant  gas  as 


it  flows  through  the  nozzle. 
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GENERAL  THLIil  ODYNAJ-'aC  EQUATIONS 


Several  basic  equations  my  be  derived  fron  a  them odynamic 
consideration  of  the  propellant  gas  mixture.  For  the  association 
reaction 


2H  *  H2  , 


the  equilibriur.  constant  is  given  by  the  expression 


Kp  ■  NH./Vr  ' 


(1) 


where  Nf  and  Nu  are  the  role  fractions  of  H  and  H,  respectively, 
•  i  ^  H 

and  P  is  the  total  pressure  in  atmospheres.  The  role  fractions 
satisfy  the  obvious  relation 


"h 


2 


♦ 


N, 


I! 


=  1  . 


(2) 


The  molecular  weight  of  the  gas  nixtur*  is  given  by  the  relationship 


K  =  £  N^.  (i  =  1,  2)  ,  (3) 

where  and  are  the  role  fraction  and  the  molecular  weight,  respec¬ 
tively,  of  component  i. 

The  expansion  of  the  propell  mt  gas  is  cassuined  to  be  isentropic. 
The  specific  entropy  of  the  gas  mixture  is  giv*n  by  the  expression 


s  -  -yr  -  P.£N.  in  N.  -  R  In  P)  ,  (;J 

where  s  is  the  entropy  in  calories  per  degree  per  gr.arr,,  K  i3  tne  mol^cu- 
lar  weignt  of  tne  gas  misture,  S°  is  the  entropy  per  ole  of  component 
i  at  one  atmospnere  and  at  a  given  temperature,  ia  the  mole  fraction 
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of  component  i,  R  ia  the  gas  constant ,  and  P  is  the  given  oressure 
in  atmospheres. 

The  specific  enthalpy  of  the  gas  mixture  at  any  point  in  the  ex¬ 
pansion  process  is  given  by  the  expression 

h*T[Z«i<»V  ^98.16.  i>]  '  (5) 

where  h  is  the  enthalpy  in  calories  per  gram,  K  is  the  molecular 

T 

weight  of  the  gas  mixture,  and  ,  and  ^  ^  are  the 

mole  fraction,  the  heat  of  formation,  and  tr.e  sencible  heat  content, 
respectively,  of  component  i. 

By  equating  the  change  in  kinetic  energy  of  the  expanding  gas  to 
its  change  in  enthalpy,  the  velocity  at  any  point  in  the  expansion  is 
found  to  be 

v  *  [v'c  ♦  2J(hc  -  h)j  1//2  ,  (6) 

where  v  is  the  velocity  in  centimeters  per  second,  h  is  the  specific 
enthalpy  in  calories  per  gram,  J  is  the  mechanical  equivalent  of  heat, 
and  the  subscript  c  refers  to  tr.e  charter  conditions. 

.  pecific  impulse  and  velocity  are  related  by  means  of  ti.e  simple 
expression 

I  ■  v/g  ,  *  (7) 

where  I  is  the  specific  impulse  in  seconds,  v  is  ti.e  velocity  in 
centimeters  per  second,  and  g  is  the  acceleration  of  gravity.  This 
concept  of  specific  impulse  assumes  complete  or  perfect  expansion  of 
the  propellant  gas. 
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MASS  FLOW  EQUATIONS 


A  fundamental  relationship  between  nozzle  configuration  and  pro¬ 
pellant  gas  .mass  flow  nay  be  derived  from  tue  following  considerations. 
The  time  rate  of  mass  flow  p  in  grams  per  second  nay  be  defined  by  the 
expression 

H  -  pAv,  (8) 


where  p  is  toe  densit  of  tr.e  gas  in  grams  per  cubic  centimeter,  A  is 
tne  cross-sectional  area  of  the  nozzle  in  square  centimeters,  and  v  is 
the  velocity  in  centimeters  per  second.  The  ideal  gas  equation  may  be 
written  in  the  form 

P  •  -ft-  KT  ,  (91 


where  P  is  the  absolute  pressure  in  atr, 03pheres,  p  is  the  density  of 
the  gas  in  grans  per  cubic  centimeter,  M  is  U  e  molecular  weight  of 
the  gas,  R  is  the  ga3  constant,  and  T  Is  the  absolute  temperature  in 
degrees  Kelvin.  Equations  (8>  and  (/)  may  be  combined  to  give  the  ex¬ 


pression 


pHT 

TfG  ' 


(10) 


for  th-  area  associated  with  a  giv«n  mass  flow  rate  at  any  point  in 
the  nozzle.  The  correspond  in g  cross-sectional  radius  is  g_ven  by  the 


express: on 


r  = 


nHT 


1/2 


(ID 


From  the  equation 


2  2  2 
x  4  y  -  rc  < 


(is; 
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which  defines  the  spherical  section  of  the  nozzle  depicted  in  Fig.  1, 
and  the  fact  that  y  =•  r,  the  following  expression  is  obtained  for  the 
distance  downstream  fror  t  e  nozzle  entrance  thn  ugh  the  spherical 
section 

2  2  1/2 

X  -  (rc  -  r  }  .  (13) 

An  expression  for  the  distance  downstream  through  the  rounded 
throat  section  :aay  be  derived  fror.  the  e  pat  ion  w.  ich  defines  the  shape 
of  the  throat.  In  terr  s  of  ar.  x,  y  coordinate  system  having  its  origin 
at  ti.e  cent-r  of  th.e  nozzle  -ntrance^  the  equation  of  sue!  a  defining 
circle  is  sir-ly 

(x  -  xt)2  ♦  (y  -  yt)2  =  ?r^  ,  (la) 


where  and  y^  are  the  coordinates  of  the  center  of  the  circle  whose 

1/2 

r  ldius  is  2r.  or  2  r  .  Reference  to  Fig.  1  shows  that 
u  C 


y, 


V 


,1/2 


and  that 


*i  *  V 


t^/< 


(13) 


(1M 


where  is  the  abscissa  of  the  point  at  station  1  where  the  circles 
defined  by  Eqs.  (12)  and  (1.*)  ire  tangent.  By  combining  Kqs.  (I4;  .nd 
(13/  and  remembering  that  y  **  r,  the  following  exnression  for  the  dis¬ 
tance  downstream  through  the  throat  section  is  obtained: 


-  fe.2  /  ..  /-1/2.2  1  1/2 

lt  *  [lTc  (r  -  !rc'2  J  J 


a  n 


In  Eq.  (1/)  tie  •  inus  sign  is  valid  in  t  ■*  convergent  s-ction  of  the 
nozzle  ,ar.d  the  positive  sign  in  tie  diverger. t. 

1  spection  of  Fig.  1  shows  tnat  the  expression  for  th.e  distance 
dovnstr  -ai  fr  >  station  2  io  g  iv-n  by  tie  relation 


9  - 


x  =  (r  -  r2)  cot  0  ,  (18) 

wh*re  r^  is  t,h»  radius  of  tr e  throat  at  station  2  where  the  divergent 
cone  of  half  angle  6  is  tangent  to  the  circle  which  defines  the  shape 
of  the  U  roat. 

Further  inspection  of  Fig.  1  indicates  the  following  particular 
relations.  The  angle  a/  associated  with  the  point  of  tangencv  at 
station  1>  is  giv»n  by  equation 


cos  a 

so  that  a  °  23.5 16°. 

Thus  at  station  1 


1/2) 

a1/  *• 


) 


0.87368  , 


and 


r  cos  a  =  0.37868  r 

c  c 


x^  *  r  sin  a  *  C.477<*1  r„  . 

1  c  c 

At  the  throat 

/  l/2 

r,  -  r  2  1  -  0.70711  r„ 

t  c  c 

and 

x.  *  x_  ♦  2^/^r  sin  a  -  1.1526  r  . 
tic  c 

At  station  2 


(19) 


(20) 


(21) 


(22) 

(23) 


r2  -  (3  -  2  cos  13°)  rt  -  0.75529  rc  (24) 

and 

x_  »  x.  ♦  2^'  “  sin  15°  r  *»  1.5186  r  . 

O  C  C 


(25) 
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KINETIC  E  QUATIONS 


The  recombination  of  hydrogen  -it  ims  to  f orr  hydrogen  molecules 
:ray  hr  represented  1  y  the  triple-body  collision  reaction  equation 


Z  ♦  2H 


H.  ♦  Z  , 


(2b) 


wher*  kj.  nc  kv  ir»  t  -  forward  and  backward  reaction  rate  constants, 
pectiv*ly ,  and  Z  is  in;,  third  body*.  if  jVl  denotes  the  concentra 
lion  of  H  in  >le-  per  cul  *c  eer  1 1  me t -r ,  t }  en  the  ti">e  rat*  of  cn..n;  e 
of  f  r  a  process  in  w;  •  t.  •  voluj  e  V  is  c.  anging  is  <.  iven  by 

t  ,  *  expre-  ion 

i  fH] 


[Hi  _£V_ 

dt 


( 


i 


dt  “  kr  LH1  M  -  hfH;]W  •  A 

Ihe  validity  of  th*  third  *  er;  r.  the  right-  » and  side  of  tv;.  (27 j  ir<y 
v-e  e  ,s .  ly  ver.fi.ed  by  a  c  uisirieratior.  of  t*  e  identity  [lij  «  ra^/V, 
wr  ere  is  the  na'.b»r  of  :  oi*s  of  i!  on  the  volume  V  of  expanding  gas 
in  t . , i s  ^tody  it  s  ore  c  nv»r.ier.t  to  d *•  a  1  vith  the  concent  of 
role  fraction  N  t  an  with  that  of  solar  concent-at ion  [li]  .  by  rak¬ 
ing  ne  r  tr.e  id *r titles 

L’J  -  NjjP/HT  (2.2) 


■  nd 


10  ■  F/UT  , 

.  (27)  ra  1-  written  in  to*  forr 


(29) 


dN, 


Xppp/tlT)-’  -  k,NH  (I  •(?)] 


♦  N„  -TT-  1"  O'V/HT  . 


♦Co.- pare  this  treat.'.e-  ♦  t  a. at  of  P-nser  . 


(30y 
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For  a  i:as:  of  r  pra'ac  of  ,;«r  the  id-al  gas  law  nay  be  written  in 


th*  fora. 


o  U  it 


PV  =  rHT/K  , 


d  n-. 


in  (IV/T)  =  -i-  In  (rrJtA*.) 


1  dK 

t  nr  • 


Hence,  t,g.  (  jC )  !  <*c ores 


TT“ 


![kfN^(P/UT):2  -  kbNn  (F/l-.Tj]  - 


H  IK 

T“  — 


My  using  h:s.  (2)  and  ())  .an  '  t‘ e  value  1.0o8  for  t;.e  atonic 
we;  -Ft  of  hydrogen,  the  expression  for  t:.-  jjecul.ar  weight  of  the  gas 


ixture  becor'»*3 


y  =  1.008  (2  -  nh>  . 


u;u  t  'on  (  J3)  thus  reduces  to  th*  following,: 


dN 

nr 


-  (2  -  NhJ  [kfN^(P/UT>r  -  kbNn  (P/mTj  . 


Fince  tf.e  e^uilibriun  ccn3t  a  t  K  defined  by  h  j,  (1)  nr. :  ti.e  **  ite 
Const;.  ;ts  in  -  O’o)  are  relit'*;!  by  t:.e  ^xur^sc  ion 


K  iiT  =  k./k,  , 

D  f  b  ' 


c.  !.  (  !3 i  fr  l-  written  ir.  t:  r  forr. 


k&'  -  V 

(KTV 


c«;,p'  -  WH  rA  >  . 


Finally,  since  tl  e  v-Ioc  '/  v  at  any  lis'arce  d.n.nstr*  is 


given  oy  the  relatior 


<  i  A.  /  G  L/  j 


tf.e  rate  of-  change  of  wit'  distance  i*-  given  )  y  t.  •  expression 

rt'n  kf r  -  V 

-  sir  •  -  %  *V  • 

‘  i  a  .  ;  2 


(39) 
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RADIO  DATA 

The  pertinent  th**r  ©dynamic  data  u?*d  in  thin  study  were  taken 
from  the  National  nu**eau  of  Dtandards^; .  T;,ey  include  the  total 
enti  alpy  above  the  reference  temperature  293.16°K  and  the  absolute 
entropy  for  the  cl  enical  species  H  and  H„  and  t  e  logaritiurs  of  the 
equilibrium  constants  for  trie  reaction  2!!  *  H  at  te  .Der-itures  from 
160C°K  to  J>CO°K.  The  heats  of  formation  for  H  and  H  at  h9B.l6°K 
are  52,069  and  0  cal  mol  \  respectively. 

The  pertinent  oiiysical  constants  i.ive  tie  following  values: 
the  acceleration  of  gravity  g  =  960. 065  cr.  sec  '  ;  the  gas  constant 
R  -  32.0567  cm^  atm  deg  ^  rol”^;  the  i..ec  .ar.ical  equivalent  of  heat 
J  =  4.16*0  abo  joules  cal  ^ . 

The  forward  reactior  rite  constant  V.  ^  in  Kq.  (26)  is  assumed 
to  have  the  value  10^  crr,°  rol-^  sec”^  over  the  entire  range  of  tem¬ 
peratures  tnd  pressures  covered  by  this  investigation.  T  .is  assump¬ 
tion  is  made  (a)  becaus-  tie  chemical  literature  is  singularly  lacking 
in  Kin'-t.c  dita  for  the  hydrogen  recoi  blnat  ion  reaction*  and  (b)  be¬ 
cause  the  primary  purpose  of  tiiis  study  is  to  si.ow  l.ow  a  rate  equation 
1  ay  be  <iealt  wit.:,  in  a  propellant  gas  expansion  problem. 


In  tl  is  connection.,  see  Refs.  5  and  6. 


-  13  - 


COMPUTATIONAL  PKOCtDIJRK 

In  general,  the  composition  of  the  expanding  gas  may  be  deterrined 
at  any  te:  perature  and  ;  res  sure  by  the  simultaneous  solution  of  L\3.  (1) 
and  (C).  The  er  tropy  of  the  system,  the  expansion  of  which  is  assumed 
to  be  isentropic,  may  1  e  d**teri  in~d  fr>r  tie  initial  Coanber  conditions 
ty  n*&ns  of  rig.  (,+  ).  Once  the  value  of  the  entropy  l  as  b-*en  established, 
ho.  (  *)  "ay  be  ised  to  deterrine  the  pressure  at  any  temperature  ind  con- 
position.  T|,e  velocity  of  t/.e  expanding  gas  at  any  noint  r.ay  be  computed 
by  means  of  ii.  (6).  The  chamber  velocity  ;:ia-  be  determined  fror  tne 
following  considerations. 

If,  for  a  given  as:  flow  rat*  a  plot  of  area  A  versus  tempera¬ 
ture  T  is  made  for  various  te  peratures  hetwee:  T  and  T  ,  a  smooth  curve 

c  c 

having  a  well-defined  minimum  is  obtained.  This  po^nt  of  rninimun  irea 
corresponds  to  !be  t.hroat  of  the  nozzle.  lie  chamber  velocity  may,  tr.ere- 
for*,  be  confuted  fror.  E  (10;  whe-i  the  uesign  condition 

Ac  3  2  At  (aOi 

is  satisfxed. 

Jince  A.  can- ot  he  cr  put*d  "x-i  icitlv.  a  tentative  v  Ti  .e  for  A 

t  -  1  c 

.ust  he  33  u-ied  on  tre  b  is.s  •'  some  •  on::ero  cr  .i  l*r  velocity,  i  hen 
values  of  A  m.a"  be  commuted  for  virions  values  of  T.  r'ro:  a  ulot  of  these 
data  a  t.nroal  .a  re  a  corre  por.d  i  •  g  to  *  h.»  ismu-i  r  inter  ar-'a.  :auy  be  deter¬ 
mined.  If  t...  i.rocesr  .5  r*p*at.e  i  for  several  val  *s  of  A  nd  if  a  plot 

of  A  /A.  v-rsjs  A  is  made,  hi, eg  Vv  interpolation.  val  es  of  A  nd  A. 

c  t  c  '  J  c  t 

nv  b-  ohta  red  which  sat  s  f  y  h  ..  (^0 1.  This  r<  c-dure  far  (determining 
A  ,  A^,  a;  a  vc  is  perfectly  g-ne  mil  ar.a  ,3  independent  of  t  .»  type  of 


flov  considered 


The  ie.ngt.  of  the  nozzle  ::ay  be  confuted  fro;;,  a  knowledge  of  end 
conditions  only  by  rre.an.-;  of  jIjS.  (13),  (17),  and  (Id)  for  constant  com¬ 
position  flow  and  for  instantaneous  e  1  cal  e^uilibr-’u..  flow.  The 
kinetic  case,  however,  requires  a  * '-pause  roced  .re  which  involves  the 
integration  of  rM.  (39)  a.  well  («>  nr  e  use  of  r.  1s.  (13),  (17),  and  (H, 
For  this  purpose  . t  is  :orvenient  to  use  * r.  perutnre  as  a  narareter  and 
to  perforr  tue  intepriticn  over  te;  peratur*  intervals  wi  ich  are  snail 
enoupr,  to  ens  .re  a  depr*e  of  accuracy  c  r  3irter;t  wit;,  the  data  us»d. 

Tie  two  rep  us  where  art'Cular  car*  rust  e/.ercised  are  (a)  at  the 
v-r.  e- trance  to  t  e  nozzle,  where  f  e  ♦  rr  ;  <*r;it  ;re  gradient  is  1  r axi¬ 
n’  ,  cud  (b;  it  let  roat,  v.”  er*»  t :  ~  para  eters  rurt  k#  det-r  ir.ed  ;  re 
' .  :»*la  In  ord"r  to  s  t.sf;  tie  d-’sipn  c  »nditi  ,ns, 

.  r.-  i  atepruti  or  f  u  ( 39;  ?•  a;,  1  <•  accomplished  by  vritir.p  _t  in 
4  ■ *  form 

rtrt,/r:L 

ilVdT  *  >  (-lj 


w  ere  L  denotes  the  ler.  id.  of  the  r.ozzle,  and  :y  :u  j**ctir:y  it.  to  an 
i  ter  it .  v*  process  whirr  •  se  t!  e  net!  od  or.O  '  just  ions  the 

n  .r-r.cal  c  llc  .lu,  for  ey.nl  intervals  of  tr.e  nrpn'  ert  T.  A  list  of 
such  j.atior  .  .  5  gi  ver;  ty  ?'  1  •  n»  ^ y .  For  t*  •  present  *  e  following; 
four-po'nt  int»pratior  Torn*  .  n  s**t '  sf  tCtory : 


N  ♦  ( 


3N,' 


3N' 


U2) 


In  t ;  is  e  puati  n  ;  3  tie  t  -'..per  tu-e  •  nt»rv  1  n.  t/.e  pr'  •  d  quanti¬ 

ties  ar*  t’  e  first  derivatives  wit4  resn°ct  to  te:.ap*rat  ir-  of  the  func¬ 


tion  N 


'  fo  .r  consee.tiv"  ter  ner.atures. 


In  order  to  c.rry  out  'he  integration  '  roces'--  effectively,  it  1 3 
expedient  to  construct  tables  of  dh^/aL  and  L  values  for  a  judicious 
spread  of  .'{  vilues  at  tr»  various  ter  peraturer  considered,  starting 
’.aloes  of  N  may  be  obtained  by  reference  to  those  c  Iculaterf  for  the 
instantaneous  e  juilit  riur  case. 

Once  a  set  f  values  of  N  is  obtained  for  the  temperature  range 

H 

frorr.  T  to  >:re  arhitrur.'  ?  .  correspor  d an  values  of  dN,,/dL  and  of  L 
c  e'  '  H 

ray  V  e  obtained  fro’  tie  ;  reviously  constructed  tables  by  interpolation. 

Values  of  dl./dT  r  ay  e  derived  fro.:  ,a  plot  of  L  versus  T  or  by  means  of 

an  ary  ropriate  differentiation  formula  as  given  in  Ref,  7.  ^,’Jiti.on  (*2) 

r  i.v  then  he  used  to  correct  ire  original  rough  values  of  N„  as  often  as 

ri 

necessary  to  achieve  convergence. 
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RfcibULiT  ANb  DLX'I^IOri 

The  ’  ole  friction  composition  of  the  gas  entering  the  nozzle  at 
r  temperature  of  3500°K  in  a  pressure  of  20  atr.  is  0.1235  H  and 
0.8?6>  H p,  and  its  specific  entropy  is  22.918  cnl  deg  ^  gro  \  with 
these  data  as  initial  conditions,  a  mas  flow  rate  of  1000  gm  -^ec  \ 
the  nozzle  configuration  of  Fig.  1,  nr.d  a  2-to-l  chair.ber-to-thro.it- 
area  ri  tio,  U  e  adhat  at  c  flow  of  hydrogen  to  an  exi.aust  pressure  of 
1  atr.  w,ts  studied,  assuring  firr.t,  constant  cor  position  flow,  then  in¬ 
stantaneous  chemical  equilibrium  flow,  ami  finally,  kinetic  chemical 
equilibrium  flow.  I  he  r 'suits  are  presented  graph. c.lly  in  Fi43.2 
tiirou  h  6. 

Inspection  shows  that  the  shape  of  the  kinetic  flow  nozzle  is 
identical  with  that  of  the  instantaneous  flow  nozzle.  This  is  due  to 
the  fact  that,  since  the  magnitude  (10^)  of  the  forward  rate  constant  kf 

is  for  all  practical  purposes  infinite,  the  chamber  parameters  are  numeri¬ 
cally  equal  (as  fir  as  significant  figures  are  concerned)  for  tt.e  kinetic 
and  t’.e  instantaneous  equilibrium  flow  cases.  .  ince  the  chamber  para- 
i  »t»*rs  determine  the  geometry  of  the  nozzle,  it  follows  th  it  the  shape 
of  the  nozzles  must  be  identical. 

As  uas  t »  be  expect'd  the  kinet.c  recults  are  for  t  e  ,ost  part  in¬ 
termediate  between  the  cor. s’ ant  composition  ar  d  t  -  instantaneous  equili¬ 
brium.  flow  results.  This  fact  was  invaluable  in  estimating  starting 
values  for  the  iteration  -rocess  used  to  in' egrate  hq.  (39;.  A  tempera¬ 
ture  interval  of  50°K  was  er  ployed  throug:  out  t:  e  integration  process.  A 
smaller  .nt*rvil  might  have  ,  .ven  smoother  results,  but  t  e  increased 
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labor  would  have  been  greater  tnan  the  accuracy  of  either  the  initial 
data  or  ty  e  final  res  .Its  would  warrant. 

Figure  6  depicts  the  nozzle  configuration  for  the  three  types  of 
flow  investigated.  It  is  interesting  to  note  (a)  that  instantaneous 
equilibrium  flow  requires  a  longer  (and  larger)  nozzle  U an  does  con¬ 
stant  composition  flow  tnd  (b)  that  for  kinetic  equilibrium  flow  the 
nozzle  is  slightly  shorther  than  for  instantaneous  equilibrium  flow. 

It  may  also  be  noted  tn.it  nozzle  design  affects  specific  impulse 
in  the  following  manner.  The  specification  of  a  2-to-l  chai.ber-to- 
throat-area  ratio  results  in  a  charier  sp-cific  impulse  of  l4l.^  sec 
for  t.ne  con  tar.t  composition  flow  nozzle  and  I36.O  s**c  for  both  the  in¬ 
star. taneous  -nd  the  Kinetic  equilibrium  flow  nozzles.  It  is  obvious 
f rorr.  cq,  (10)  that  a  zer)  chamber  velocity  would  require  an  infinite 
chamber  area,  a  factor  which  is  not  fe,a  ihle  in  practice.  Furthermore, 
a  zero  ch  ur.ber  velocity  vould  yield,  at  an  exhaust  pressure  of  1  ntn,  a 
specific  impulse  of  82). 2  sec  in  te,id  of  335.-*  for  constant  conpo- 
sition  flow  and  373*3  sec  instead  of  884.3  for  instantaneous  equilibrium 
flow.  It  is  thus  evident  that  nozzle  cmfigurat  on,  specifically  chamber- 
to-thro  it-are  i  ratio,  has  a  dep:nite  effect  on  the  magnitude  of  specific 
in. pulse  or  jet  velocity. 

In  the  present  example  the  effect  of  cuemic.l  kinetics,  l.e,,  the 
rate  of  recombination  of  hydroger  atoms,  on  specific  irnulse  is  not  very 
prono jr.ced ,  as  the  figures  88.+.  )  and  '7^ . 6  at  an  exhaust  pressure  of  1 
atm.  for  instantaneous  ind  ;inetic  eviilibrur'  flow,  respectively,  indi¬ 
cate.  T  ^  i  3  is  due,  of  course^  to  the  a  an:*  ;d»  of  the  forward  reaction 
rate  constant  k^  and  to  the  nature  of  the  chemical  species.  . r.  ot.her 
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expanding  propellant  gases  in  which  the  rate  constants  associated 
.  ’t  the  pertinent  che  deal  re  »ctiot.s  are  *iu  uerically  smaller  than 
the  hydrogen  rate  constant  ir  d  ‘ .  e  ci  e. ideal  species  themselves  are 
less  active  than  hydrogen,  t- e  general  effect  won  Id  be  to  displace 
tie  kinetic  res  .Its  toward  the  constant  composition  flow  values. 

In  such  cases,  even  if  several  differential  e  juatio:  s  similar  to 
(IV)  are  involved,  t.r.e  method  of  integration  indented  in  this 
pa D'T  would  l'  valid,  but  the  computation  would  best  be  performed  on 
automatic  equipment. 
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THKOICH  A  HOCKKT  NOZZI.K 
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ROCKKT  NOZ/. IiK 


CONCLUSIONS 


7 r.e  following  conclusion:-.  ray  h<-  dr  iat.  fro;.  the  foregoing  study. 

1.  Ti.e  effect  )f  c.  .epical  kinetics  on  n-c.et  pro pell ant  g-a  s-s 
»xp  .ndiny  through  a  .uozzl*  has  t»-»n  determine. i  t  v  ti  e  *r.  jit  --meows 
fioluti  on  of  thr***  sets  of  e  oi.it  ions  base  :  on  ti.er  odynu.  *c,  nozzle  con¬ 
figuration,  .r.«i  c  eirdcnl  n*.ctl  >n  rate  consider  it,  ions. 

T.ue  fundamental  aes^pr.  c  ilcu^at:  t.hI  n-4iire--:it,  4  idf  fr. 
t.‘  e  r  i.  flow  rate,  is  t: .•  '•pec  ricition  of  a  finit*  c  a  .r -r-to-tun.  t- 
\T'  h  ratio. 

3.  Nozzles  designed  on  tne  ksi.;  of  ct  e-  deal  kinetics  for  perfect 
expar  iion  are  a  or  ewf.nt  nailer  t  an  inct  ant  arieons  e  pilibriur.  flow  noz¬ 
zles  but  ri’-e  lnrg-r  t  r,  cons'. ant  co-*  position  flow  nozzles,  t1  *•  d.ff'T- 
-•  ee  0-- enoin-'  on  t;.e  n.agnit  ide  of  l'  e  reaction  rat-  Involved. 

t\.  n.»  -ethod  of  step. due  integration  d-  crlb»d  ray  be  npplieu  to 

pr  pell  nt  uysteff,  wbic  inv'ilv*  several  differ'*:  *i  i  e  .  .ti  ns  witt  ,  of 
course,  1  correspond  lry  inc-ca  *•  lr  c  ■:  putational  labor. 
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